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ABSTRACT 

Conjugate  basis  functions  which  are  continuous  through- 
out a  model  are  used  to  approximate  a  stress  field.   This 
conjugate  approximation  has  less  mean  error  than  approxima- 
tions calculated  using  conventional  finite  element  methods 
and  is  a  better  approximation  at  extreme  values  of  stress. 
This  thesis  consists  of  two  major  parts.   The  appropriate 
analytic  expressions  for  conjugate  stress  fields  are  set 
forth  and  applied  to  the  constant  strain  triangle  finite 
element.   Analyses  are  then  performed  with  these  new  stress 
calculations  and  the  results  compared  to  some  previous 
stress  calculations. 
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I.   INTRODUCTION 

A.   BACKGROUND 

In  a  recent  paper  Oden  and  Brauchli  [Ref .  1]  state  that 
one  of  the  difficulties  encountered  in  using  the  finite 
element  method  to  approximate  the  stress  field  in  a  model 
has  been  that  the  stresses  calculated  are  at  best  only 
averages  over  each  element  and,  in  general,  are  discontinu- 
ous between  elements.   Continuous  stresses  can  be  obtained 
only  if  strain  degrees  of  freedom  (DOF)  are  introduced  at 
each  node.   This  is  achieved  at  great  computational  expense. 
Some  of  the  schemes  which  are  used  to  calculate  the  state 
of  stress  at  a  node  are  mentioned  in  Ref.  1  and  another 
method  which  is  based  on  conjugate  approximations  is  pro- 
posed.  In  Ref.  2  Oden  expands  on  this  method  explaining 
the  theory  of  modeling  a  function  in  general  and  using  con- 
jugate space  functions  for  modeling  in  particular. 

Briefly  the  methods  of  evaluating  nodal  stresses  re- 
ferred to  in  Ref.  1  include  taking  an  average  of  the 
stresses  of  each  element  incident  at  the  node;  using  a 
weighted  average  method  developed  by  Wilson  [Ref.  3];  using 
averages  based  on  equating  nodal  forces  and  element  stresses 
in  a  procedure  developed  by  Turner,  Martin  and  Weikel 
[Ref.  4];  or  using  stiffness  matrices  proposed  by  Gallagher 
[Ref.  5]  which  are  developed  relating  both  the  stress  and 
displacement  fields  of  each  element.   Among  all  methods,  the 
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proposed  method,  using  conjugate  bases,  has  been  shown  to 
give  the  best  fit  in  a  least  squares  sense. 

B.   SCOPE  AND  OBJECTIVES 

In  this  thesis  the  conjugate  approximation  theory  is  set 
forth  in  detail  and  is  applied  to  illustrative  examples. 
The  appropriate  analytic  expressions  for  conjugate  stress 
fields  are  derived  for  the  constant  strain  tiangle  (CST) 
finite  element  and  incorporated  into  a  CST  computer  program. 

Numerous  examples  of  in-plane  problems  are  worked  com- 
paring both  the  CST  solution  and  the  conjugate  basis  method 
(CBM)  solutions  and  the  exact  solution  if  it  is  available. 

Finally,  some  conclusions  are  reached  concerning  the 
type  of  in-plane  problems  where  the  CBM  may  be  used  success- 
fully and  what  the  cost  in  additional  computer  time  may  be 
for  this  continuous  best  solution. 
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II.   CONJUGATE  APPROXIMATION  THEORY 

A.   MODELING  A  FUNCTION 

A  function  F(x)  over  a  region  or  domain  R  can  be 

modeled  by  another  function  F(x)  in  the  domain  R  where 

F  ~  F  and  R  ~  R.   The  region  R  is  subdivided  into  elements 

where  each  element  has  two  or  more  nodes  depending  on  the 

degree  and  type  of  the  approximation  function  over  the 

element.   Within  R  all  elements  have  the  same  number  of 

nodes  and  each  node  in  the  global  model  R  is  designated 

X  ,  (A  =  1,G),  where  G  is  the  total  number  of  system  nodes. 

The  domain  of  the  element  e  is  r   and  each  node  within  this 

e 

N 
local  model  is  designated  x  ,  (N=l,  N  ).  where  N   is  the 

°  e        '   e  e 

number  of  nodes  in  an  element. 

For  example,  a  one-dimensional  axial  bar  may  have  the 
following  three  element  representation   (here  R  =  R) . 


- 

0  " 

L, 

12     3     4 
X     X     XJ    X 


1     2 
x     x 

e     e 


r-,    r„    r  1  r  "* 


Figure  1.   Modeling  a  Bar. 


The  value  of  F(x)  at  X   is  denoted  F  ,  that  is 
FA  -  F  (XA). 
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Over  each  element  the  function  F(x)  is  represented  by  a 

local  function,  f  (x) ,  which  is  only  defined  on  the  element 

N  N  N       N 

e.   Its  value  at  x   is  denoted  f  ,  that  is,  f   =  f (x  ) . 

e  e  e       e 

Hereafter,  the  index  e  denotes  element  e.   In  the  figure 

N 
below,  the  functions  F(x),  and  f«  for  element  2  (i.e.  e=2) 

are  shown. 


F(x) 


X. 


Figure  2.   Global  and  Local  Functions. 

B.   INTERPOLATION  FUNCTIONS 

A  set  of  independent  local  interpolation  functions 

e 

ipN(x),  (N=l,  N  )  ,  is  defined  over  each  element  such  that 

ipN(x)  =  1  at  node  N  and  zero  at  all  other  nodes.   (For  ele- 
ments with  linear  displacement  fields,  such  as  the  elements 

used  in  this  thesis,  \b„    is  linear  between  nodes.)   For  a 

TN 

one-dimensional  2  node  element,  a  set  of  interpolation  func- 


tions for  representation  of  a  linear  field  is, 

♦J  -  i  -  C  ;  ^2  =  c 


(i) 


where  £  is  a  local  coordinate  over  the  closed  domain  [0,1] 
and  origin  at  local  node  point  1. 
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Figure  3.   Local  Interpolation  Functions 


The  local  function,  f  (x) ,  is  approximated  over  the  ele- 


re 


ment  by  f  (x)  where 


fe(x)  =  fl    l[>*(x),  N  =  l,  Ne. 


(2) 


—  e         N   e 
[The  convention  f  (x)  =  f   i^N(x)  means  summation  over 

repeated  indices,  that  is 


re 


f^-(x)  =  fjK(x)  +  ff^(x)  +  ...+  f*Je  *JJ  Cx). 
er  1        e  2  e    N 


This  summation  convention  will  be  used  throughout.] 

—  e         N 
[Note  that  at  a  node  the  value  of  f  (x)  is  f  . ]   Here 

2  -2 

the  local  function  f  (x)  is  modeled  by  f  (x)  with  linear 

2       2 
interpolation  functions,  ty      and  i|^_  . 
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fn 

e 

n+1 


(x) 


(x) 


n+1 


-x 


Figure  4.   Approximation  to  Local  Function 


C.   INCIDENCE  OPERATORS  AND  ASSEMBLY 

eA      eN 
The  function  A.T  and  Q.  are  known  as  incidence  operators 

N       A 

and  are  used  in  relating  the  global  model  R  to  local  ele- 
ments r  .   They  have  the  value  of  1  or  0  as  follows: 

e       J 


1   if  node  A  in  R  is  coincident  with 

node  N  in  r   in  the  assembled  model, 
e 

0   otherwise 


(3a) 


1   if  node  N  of  element  r   is  coincident 

p 

eM  with  node  A  of  the  assembled  model.        (3b) 

fiA   = 

0   otherwise 

A  is  called  a  mapping  function  because  it  gives  the  value 

of  a  local  node  to  the  incident  global  node.   Likewise,  0, 

is  called  a  decomposition  function  because  it  gives  the 

global  value  to  the  corresponding  local  node. 

F   =  A.T  f   (Note  this  is  summed  over  both  (4) 

N    e  J   M   A 

e  and  N . ) 

,N    ^  ^A  /cN 

fe  =  fiA  F  (5) 
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The  model  of  a  function  F(x)  in  terms  of  jJj   element  in- 


terpolation function  can  now  be  given  by 
F(x)  ~  F(x)  =   Y^  ^U>  "ApA 


e  =  l 


E=l,  N 

: 

N  =  l,  N 

< 

A=l,  G 


(6) 


where  N   is  the  number  of  elements. 

r 

Equation  (6)  is  true  everywhere  except  at  the  global 

nodes  with  m  elements  coincident.   Here  the  value  of  F(x) 

is  mF(x)  times  greater  than  it  should  be.   Accordingly,  at 

the  nodes  Equation  (6)  becomes 

E 


(7) 


e  =  l 
We  note  that  Equation  (7)  gives  continuity  of  F  at 

nodal  point s . 


A  set  of  global  interpolation  functions  4>.(x)  are  now 


defined  as  follows: 

4>A(x) 


E*n 


(x)   n 


N 


/e^l 


1 
m 


£ 

e  =  l 


N 


♦  „<*>  fiA 


except  at  nodes 


at  nodes         (8) 


It  will  be  demonstrated  in  the  examples  that  <j>A(x)  is 
a  linear  combination  of  the  element  interpolation  functions 
ip  '  s  .   Moreover,  since  the  local  i|i   are  independent,  the 
global  (J),  functions  are  independent  also.   There  is  a  unique 
(J),  for  each  global  node.   If  the  element  interpolation 
functions  are  independent,  the  $A's  are  independent  func- 


tions .   Now 


F(x)  =  <J>A(x)  FA  =  ^2      ^(x)  f* 

e  =  l 


'N 


(9) 
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The  domain  of  a  global  interpolation  function  <j>.  extends 
over  all  elements  incident  at  X  .   This  function  has  the 
value  1  at  X   and  zero  at  all  other  nodal  points.   For  ex- 
ample, with  the  previous  linear  interpolation  functions, 
and  the  three  element  axial  domain,  we  have 


1  - 


**"*  X 


Figure  5.   A  Global  Interpolation  Function. 

D.   FUNDAMENTAL  MATRIX 

The  set  of  G  linearly  independent  functions  {<j>A)  form 

the  basis  for  a  G-dimensional  space  $  where  G  is  the  total 
number  of  global  nodes.   The  fundamental  matrix  of  the  space 
$  is  defined  as 

CAr    E    <cj>A,4>r>                  A,T    =    1,G  (10) 

The    notation    <<$>Ai$T>    is    used    to    denote    the    inner  pro- 


duct of  <}>.  and  (f>r.   Here  the  inner  product  is  taken  as 

<<f>A,<i>r>  =  /R  4>A  <i>r  dR     A,r  =  i,g  (li) 
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Since  /   <f>.  (J).  dR  >  0,  C.p  is  a  positive,  definite 
matrix  and  its  inverse  can  always  be  formed.   It  is  defined 


as 


(c^r1  -  c 


AT 


(12) 


This  will  be  the  transformation  matrix  used  to  form  a 
new  set  of  basis  functions  in  the  G-dimens ional  space.   More- 
over, since  /   <J>.  (J)„  dR  =  /   (J)„  <p .     dR,  C.p  is  also  a  symmet- 
ric matrix. 

It  is  convenient  to  construct  C.p  from  local  interpola- 
tion functions  as  follows.   Substitute  Equation  (8)  in 
Equation  (10)  to  obtain 

E     E 
CAF  '  <W  =I]  T!   «A  ^f   <*2-  *m" 


(13) 


e=l  f=l 


e    f 
Here  <ip.T,  ij;  >  is  the  inner  product  of  the  local  interpola- 

N    M 

tion  functions.   It  is  non  zero  only  if  e=f  since  each  func- 
tion is  only  defined  on  its  own  element.   The  inner  product 

e    e 

<}Pni    V    >    is  called  the  local  fundamental  matrix  for  element 
N    m 
e 

e»  CNM- 

Equation  (13)  now  becomes 

-  E 


'AT 


E 

e  =  l 


6N  eM   6 

"a  nr     CNM 


(14) 


E.   CONJUGATE  SPACES 

A  conjugate  basis  is  one  which  is  reciprocally  related 
and  interchangeable  as  to  properties  with  the  original 
basis.   A  set  of  G  linearly  independent  functions  on  R, 
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{(f)  .},  a  so-called  conjugate  basis  set,  must  satisfy  the  re- 


ciprocity relation, 


..A    ,    s         .A        ,  .A   _    ,1    A=T 

<<J>    ,  <J>r>    =    5p       where    Qj,    =    {Q    A^r 


(15) 


This  means  that  the  sets  {<{>*}  and  { <j)  }  form  a  countable 
biorthogonal  basis  of  the  space  $,   Any  sets  (<j)  }  and  (^a) 
which  satisfy  Equation  (15)  are  referred  to  as  conjugate 
basis  of  $ . 

In  order  to  define  { $  }  uniquely,  they  are  formed  as 
follows 

4>A(x)  =  CAr  «J)r(x)  (16) 

It  can  be  shown  that  the  $  's  are  independent  functions, 

and  that 

*A(x)  =  CAr  /(x)  (17) 

A  AK  AK 

Moreover         <$    ,  cp  „  >    =    <C       ^K><^r>    =    C         <<^K'^r> 

-  C    C  r  =  <$r   as  required  by  Equation  (15). 
Using  the  incident  operators  defined  in  Equation  (4), 


the  local  conjugate  basis  interpolation  functions  are 


(18) 


and  making  the  substitutions  from  Equations  (16)  and  (9) 

gives 

xr  e  m         AT 

(19) 


,N         "N    -AT     , 

ie  =  "A  c       *r 


*N    _AT    !?       ,f,         *M 

=  n.   c      z     ^M(x)  fir 
f=l 


A,r  =  i,g 

M    =    1,N       (20) 


Equation  (20)  highlights  the  most  important  feature  of  the 
local  conjugate  basis  functions.   The  functions  for  each 
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node  of  each  element  are  linear  combinations  of  the  basis 
functions  of  all  G  nodes.   This  gives  each  local  node  a 
basis  which  is  continuous  over  the  entire  model.   The  first 
example  gives  a  numerical  illustration  of  this. 

F.   APPLICATION  OF  THEORY  TO  FINITE  ELEMENT 

Recall    Equation    (10)    where    F(x)    =    FA(j)    =    Z       i/>®    f^     .       An 

6=1    r6 

alternate  representation  of  F(x)  in  terms  of  the  cj>   basis 
is  obtained  as  follows 


"    M   =  \    ie  fN 
e  =  l 


K 


(21) 


Here  Fv    is  the  value  of  F(x)  with  respect  to  the  con- 

jugate  basis.   It  is  calculated  by  forming  the  inner  product 

of  F  and  $     .   From  Equation  (21)  we  have 
K 

<F,V  -  <F^N,V  =  FN<^N,V 


=  F   6N 
N  \ 


=  F 


K 


(22) 


The  local  values  of  F  are  obtained  by  the  usual  decom- 
posit  ion , 


e   !a 


;a 


fN  "  «NFA  "  °H  ^'V  =  <F'V 


fN  -  at   FA 

e     N 


<F,4.A>  =  <F,^> 


(23) 
(24) 


Now  the  original  function  can  be  approximated  in  the 
conjugate  basis  using  Equation  (21).   This  is  the  procedure 
used  in  solving  problems  using  the  conjugate  basis  method 
(CBM) . 
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G.   BEST  APPROXIMATION 

Before  proceeding,  it  is  important  to  verify  that  ap- 
proximating a  function  using  the  conjugate  basis  functions 
will  give  the  best  approximation  in  a  least  square  sense. 

This  is  proven  simply  and  clearly  with  a  simple  example 

2 
The  function  F(x)  =  x   is  to  be  approximated  in  the  in- 
terval x=0  to  x=l  by  a  linear  polynomial  F(x).   Here  bases 
of  F  are  functions  $  =1  and  $9=x 

F  =  a^  +  a2<}>2  =  a^l)  +  a2(x)      (25) 

The  difference  between  the  actual  function  and  its 
approximation  is  called  the  error,  E 

E  =  | F  -  F |  (26) 

The  function  J  =  J(a.)  is  a  measure  of  the  error  over 

x 

the  interval. 

J(a.)  =  /   E   dx  =  /  [(a,  +  a~  x)  -  x  1  dx 
1     o  o     !     2 

J(a.)  -  a±2+    aia2  -  f  a±    +   j  a/    -  §  a£  +  |     (27) 

The  error  is  minimized  by  setting  the  partial  deriva- 
tion of  J(a.),  with  respect  to  a.,  equal  to  zero  and  solv- 
ing the  resulting  set  of  equations  for  a.'s. 

3J(ax)  2 

=  0  =  2a1  +  a0  -  — 

3a1  12    3 

3J(a  ) 

^—   -    0  -  a,  +  fa2  -  £  (28) 

Now  ax  -  -  j   ;  a2  =  1       and  F  =  -  |  +  x  (29) 
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Using  the  conjugate  basis  method  with  basis  functions  T\).t 
F  is  determined  as  follows: 

1.  Assume  the  local  basis  functions  are  ty       =  1  -  x 
and  \J/„  =  x. 

2.  Form  the  local  fundamental  matrix  c.T  =  <ty„ ,  U»  >  . 

NM     N  rM 


Here  we  obtain, 


CNM  =  hi    l]  (30) 


3.  Take  the  inverse 

NM     .  4  -2,  .... 

c    =   [_2   ^J  (31) 

N    eN   AT 

4.  Form  the  conjugate  basis  function  ip   =  £2.  C   ({). 

in  this  example  <£„  =  ip„  because  there  is  just  one  element. 

ij;1  =  C11  <J)X  +  C12  <J>2  =    4  -  6x 

2  21  22 

rp      =    C         <{>1   +    C         <j>2    =    -    2   +    6x  (32) 

e  e 

5.  Determine    a.    where    a.     =    <F,     \b .  > 

x  x  x 

r1       2a    a            1                           r1  2a                  1 

a      =    J       x    <J>,dx    =    Yyj       a2    =  x    ^2         =    Z" 
o                                                        o 

1    e  . 

6.  Approximate    F(x)    with    F    =    E    a.     lb  i  =  l,2 

,     x    Te 
e  =  l 

F    =    J2<4-6x)    +   i("2+6x) 

=    -    J  +    x  (33) 
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Figure  6.   Best  Approximation 

This  is  the  same  approximate  function  obtained  using  the 
previous  calculus  method.   Therefore,  the  best  approximation 
of  a  given  function  F(x)  in  the  subspace  <j>  is  the  function 
F(x)  whose  components  are  the  inner  product  of  F(x)  and  the 
conjugate  basis  functions.   Reference  2  contains  a  formal 
proof  that  the  CBM  is  the  best  fit  in  the  least  square  sense 
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III.   APPLICATIONS  TO  DEMONSTRATIVE  EXAMPLES 

In  all  of  the  problems  considered  in  this  thesis,  the 
displacement  must  be  in  the  domain  of  the  structure.   For 
one-dimensional  problems,  the  model  is  a  bar  displaced  along 
the  x  axis.   The  two-dimensional  model  is  a  plane  with  dis- 
placements in  either  the  x  or  y  direction  or  both. 


Figure  7.   One-  and  Two-Dimensional  Domains. 

Appendix  A  contains  a  brief  summary  of  the  stress- 
strain  and  strain-displacement  relations  for  the  in-plane 
problem. 

A.   ONE-DIMENSIONAL  PROBLEMS 

1  .   Nonlinear  Displacement  and  Linear  Properties 

To  demonstrate  certain  properties  of  conjugate  ap- 
proximation functions,  example  9.2  in  Ref .  2  will  be 
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explored  in  detail. 

a.   Statement  of  Problem 

Determine  the  stress  in  a  nonhomogeneous , 3  unit 
length  bar  based  on  approximate  displacement  fields.   The 
stress  is  given  by 

a(x)  =  k(x)  4"l20  (34) 


where 


k(x)  =  k  (1+x) 


(35) 


is  the  material  modulus  in  units  of  force/length   and  u(x) 


is  the  displacement  field.   Assume  that 

2 
u(x)  =  a[l  -  f  ] 


(36) 


where  a  is  a  small  constant  and  x  is  measured  in  units  of 
length . 

With  this  displacement  field,  the  exact  stress 
distribution  is 


2ak 


a(x)  = 


x(l+x) 


(37) 


The  model  R  consists  of  three  one-dimensional 

elements,  each  of  unit  length,  that  is  R  =  R  =  [0,3]  = 

3 

U  Q       . 

e  =  l  e 


(1)       (2) 


(3) 


x1       x2 


Figure  8.   Model  for  Problem  III  Al 
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The  local  basis  functions,  for  a  linear  field, 
corresponding  to  a  typical  element,  e,  are  the  interpolation 
functions  shown  with  Equation  (1), 

i|^(0  =  1  -  i  and    <|^(0  =  5 

b.   Global  Interpolation  Functions 

The  global  interpolation  functions  are  calcu- 


lated using  Equation  (9), 

3   e 


N 


<j>A(x)  =   E   ft^   ^N(x) 
e  =  l 


A  =  1,4 
N  =  1,2 


NOTE:   The  notation  (x)  will  be  omitted  unless  it  is  needed 
for  clarity. 

The  above  expression  gives, 

112  2  3  3 

11  2     1  12  2     2  13  2     3 


i  ^  +  *o 


+    0 
2 


+    0 
2 


+    0 
3 


+    0 
3, 


11  2    1 

<j>2    =    ^^1    +    ^2^2    +    ^2^1    +    ^2^2    +    ^2^1    +    ^2^2 


+    1    1^2    +    1    4>l    +    0 


+    0 
3 


+    0 
3 


2  2 

11  2    1  12  2    2  13  2    3 


+    0 
1 


+    0 
2 


+    1    ^^    +    1    ^    +    0 


2  3  3 

2  12  22  II  2     3 

4»4=    bJ*J  +  n^£  +  nj^  +  a^2  +  fi^i  +  ^2 


=    0        +0        +0        +0 

These  are  shown  in  Figure  9. 


+  0 


+  1  ip; 


(38) 
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X 


♦. 


♦. 


1- 


0 


X 


12  3 

Figure  9.   Set  of  Global  Interpolation  Functions 
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Recall  that  using  Equation  (8)  the  value  of  <j>A 
at  a  node  where  2  or  more  elements  are  incident  is  still  1. 
c.   Local  Fundamental  Matrix 

The  local  fundamental  matrix  (LFM)  is  formed 
using  the  local  basis  functions  as  discussed  in  Section  II-D. 


e  ^ i e        i e^ 


NM 


N'     rM 


N,M    =    1,2 
e    =    1,3 


c11    =    /^(SH^Od^    =    /1(l-C)(l-C)d^=    | 
o  o 

i12  -^/^con^u)^  =  /V-m  d^  =  \ 


l21  =  £y\un\(z)*K  =  /*(Ou-5><u-  \ 


c22    «    /^(U^tOdC    =    {^OtOd?      -    ~    (39) 


m       !      lr2  1, 
N°W     CNM  =  6[1  2] 


(40) 


Since  the  elements  are  geometrically  identical,  we  have 

12      3 

CNM    CNM    CNM* 

d.   Global  Fundamental  Matrix 

The  global  fundamental  matrix  (GFM)  is  formed 
using  Equation  (14), 


E   6N   eM   e 

cat  ~  En  "a    fir     cnm 

e  =  l 


N,M  =  1,2 

A,r  =  i,4 

E  =  3 


Since  C.p  is  a  4  x  4  symmetric  matrix,  a  total 
of  ten  separate  calculations  would  be  needed  to  determine 
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all  of  the  elements  of  the  GFM.   The  calculations  for  c 


11 


and  c?9  are  shown  below. 

11  11  11  11 

11    ]  12  1  2    1  1  2    2    1 

Cll    =    "l"l    dll+    "A  C12+    fiA  C21+    "l"l    C22 

22  22  22  22 

11    2  12  2  2    1  2  2    2    2 

+  njnj  c11+  ojnj  c12+  njoj  c21+  ojoj  i22 

33  33  33  33 

11    3  12  3  2    1  3  2     2    3 

+    filfil    ^11+   "Pi  ^12+    "Pi  ^21+    "l"l    ^22 

=  i[l»l-2  +0+0+0  +0+0+0+0+0    +0+0+0] 

2 
6 

11    1  12  1  2    1  1  2    2    1 

C22    =    fi2fi2    cll+    n2fi2  c12+    ^2^2  c21+    SS    °22 

22        22        22        22 

11  2      12  2      2  1  2      2  2  2 

+  fi2n2  Kl*    fl2fl2  C12+  °20i  '21+  «2fi2  c22 

+  fi2fi„  c ._+  ^2^0  c12+  ^9^2  °21+  ^2^2  C22 

=  -^[0+0+0  +2  +2  +0+0+0  +0+0+0+0] 
6 

4 

6 

The  complete  GFM  and  its  inverse  are  shown 
below . 


(41) 


2 

1 

0    0 

1 

1 

4 

1    0 

6 

0 

1 

A    1 

0 

0 

1    2 

26  -7   2  -1 

-7  14  -4  2 

=  i  rm   =  — 

'AT          0  14  1                   45      2  -4  14  -7 

-1   2  -7  26 


.AT 


(42) 
In  practice,  the  entire  set  of  calculations 
using  incidence  operators  is  avoided  and  the  LFM  elements 
are  placed  in  the  GFM  according  to  the  following  scheme. 
Each  column  and  each  row  of  the  GFM  is  asso- 
ciated with  one  and  only  one  global  node.   Likewise,  each 
column  and  each  row  of  the  LFM  is  associated  with  just  one 
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global  node.   In  assembling  the  GFM,  the  ij  element  from  the 
LFM  associated  with  local  nodes  i,j  is  added  to  the  k,l  GFM 
element,  where  global  nodes  k,l  correspond  to  local  nodes 
i,j.   If  several  elements  are  incident  at  a  global  node, 
then  the  appropriate  elements  from  each  of  the  connecting 
elements'  LFM  will  be  summed  in  the  GFM.   This  procedure 
avoids  numerous  calculations  where  the  incident  operator  is 


zero  . 


e.   Local  Conjugate  Basis  Functions 

The  local  conjugate  basis  functions  are  calcu- 


lated using  Equation  (19), 

,n,  .   J:n  at  ,  ,  . 

Jte(x)  =  fiA  C    <J>r(x) 


e  =  1,3 
N  =  1,2 

A,r  =  i,4 


The  calculation  for  node  1  of  element  1  is 


shown  below. 
1 


1   11 

1  C  ^1 


'1   c  v: 


&         =      ft^"      C"""  """(J),       +     ^      C*"4>«J      +      ^T      c~"$i      +     fiT      C~"<J> 


ij    c13cp: 


1      14 

1  c   u 


(43) 


*1       21.  Jl       22,  Jl       23.  o1       24a 

*1       31..  }.l       32.  Jl       33.  *1       34. 

+    J2_    c       9-    +    fi_    c       92    +    u       c       9„    +    U„    c       9, 

*1      41,  *1       42.  *l      43.  *1       44a 

+    fi,     c       9-1    +    "a    c       9«    +    52      c       <p      +    fi,     c       9, 

e  *1 

r-  j-ijTi  *.     i     nN  1    for    £2.. 

For    node    1    of    element    10.=  1 

0    for    all    other 

AT  1 

Therefore,  substituting  for  C    and  9A  above  $)      becomes 

±]_    =  |3-[26  <ij  -  7(^2  +  ty\)    +    2C*2  +  *J)  "  ^2]      (44) 
where  the  ^'s  in  brackets  are  the  local  basis  functions  de- 
fined in  Equation  (1)  . 
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Notice  that  only  one  row  of  the  above  calcula- 

6N 
tion  was  nonzero  -  that  is  the  row  where  fi.  =  1.   Using  this 

fact  reduces  the  actual  calculations  in  a  manner  similar  to 

the  method  used  in  forming  the  GFM. 

Consider  the  nonzero  row  of  the  above  calcula- 
tion.  Factoring  out  the  common  incidence  operator  leaves 
four  terms  which  are  the  products  of  the  Tth  element  of  the 
A  row  of  the  GFM  inverse  times  <}>„  .   Recall  that  at  a  node 
fyy    =  1  and  the  local  conjugate  basis  functions  X  =  X   be- 
come equal  to  the  A  row  of  the  GFM  inverse.   This  greatly 
simplifies  computer  programming  of  the  conjugate  basis 
method . 

Returning  to  the  example  at  hand.   Global  node  2 

is  node  2  of  element  1  and  node  1  of  element  2.   Therefore, 

X2    21  X2    21 

2    =  ^2  =  1    and  fiA  =  fiA  =  °  for  A  =  X»  3>  or  A* 

Showing  only  the  nonzero  terms  of  the  calcula- 
tion , 

t[  -  ±1  =  c21^1  +  c224»2  +  c234.3  +  c2%4 

-   fjl-7^    +    14(if>2    +    *J)-*(*2    +    ^1}    +    2^2]       (45) 
A    similar    argument    for    global    nodes    3    and    4 
gives 

±1    "   *3    "    T5[2^l    "4(,J;2    +    ^1}    +14^2    +    ^i)-7^2]  (46) 

and 

±1    =   fj[-l*J    +    2W\    +    4>[)    -Hi>l    +    *J)    +    26^]  (47) 

The  local  conjugate  basis  functions  are  shown 
graphically  in  Figure  10.   Notice  that  each  function  is 
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±{M    (jj) 


52 

4 

0 

^""^"^-^ 

3 

0 

\     1            ^ 

X1  ^"^ 

i 
2 

"""*^-r2 

-14 

±[U)=  ±\wijj) 


±22U)=   i*(x)(±5-) 


ijwciy) 


Figure  10.   Local  Conjugate  Basis  Functions 
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continuous  over  the  domain,  x,  and  that  the  effects  of  nodes 

remote  from  the  local  node  being  considered  diminish  as  the 

distance  from  that  node  increases. 

f .   Conjugate  Nodal  Components  and  Global 
Representation 

Calculate  the  conjugate  nodal  components  of  the 

element  stress  field  according  to  Equation  (23). 

o*    =  <a(x)  ,  4^(x)> 

Since  \p       are  in  terms  of  £ ,  o(x)  will  be  transformed  to 

functions  of  £  also, 


ak 


a  factor  of  - 


a(£)  =  2£  +  2£  for  0  <    x  <_  1 

=  2E,2+    6£+4  for  1  <  x  <  2 

=  2£2+10£+12  for  2  £  x  £  3 

has  been  factored  out  to  enable  compari- 


son with  the  results  published  in  Ref.  2.   Now 


o\    =    f1     (2£  +  252)(1-Dd£  =  § 
o 

o\  =  s1  us  +  2z2)a)dt  =  | 

o 

cr2  =  S1    (2?2+  6C  +4)  (l-C)d?  -  |i 
o 

a2  =  Z1  (2?2+  6?+  A)(C)d5  =  I2 
o 

o]    =  71  (2^2  +  10^+12)(l-5)d5  =  *£- 
1   o  6 

c*  =  f1    (2C2+10C+12)(C)dC  =  ~ 


(48) 


Now  the  global  stress  field  can  be  modeled  using  Equation 
(21),  3 


cOO    2  a*  ^(x) 
e  =  l 


N  =  1,2 
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At  the  global  nodes  we  have 


oCO.) 


3  52    1_   1_4    19  14    27  4     47  4 

6  15    6  15    6   15     6  15     6  15 


59  2_ 
6  15 


=  -  0.33 


„,,  .    3  (-14)    28.7    19.    8  .27  .  47.    4  ,59. 
0(1')  =  6  —IT-  +  15  (6  +  6~)  "  15  (T  +  ~Z>    +   TJ^ 

=  3.43  (49) 

i 

Similar  calculations  are  made  for  x=2  and  x=3. 
The  results  are 
0(2.)  =  11.7 
0(3.)  =  23.6 

These  results  are  plotted  on  Figure  12  and  are 
also  shown  in  Table  1.   Since  the  exact  stress  field  was 
known,  a  comparison  of  accuracy  is  made. 

g.   Conventional  Finite  Element  Solution  (C.F.E.) 
Here  we  suppose  the  continuous  displacement 
field  U(x)  is  approximated  by  the  piecewise  linear  field, 
U(x)  where 

U(x)  =  |[9-(j>1(x)  +  8<J>2(x)  +  5<J>3(x)].         (50) 

The  <J) '  s  are  shown  in  Figure  9.   This  linear 
field  was  constructed  to  give  U(x)  =  U(x)  at  the  nodal 
points  by  evaluating  U(x)  at  the  nodes.   The  strain  is  now 


dU(x) 
dx 


a 

9 


d<J>. 
dx" 

d<|>. 
dx" 

d<|>. 

dT 


+   8 


+    5 


d*. 

dT 

d<J>. 

dlT 


-    1         0    <    x    <    1 


=    3         1    <    x    <    2 


=5         2    <    x    <    3    (51) 
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The  conventional  solution  is  obtained  from 

/  n    i  /  n  dU(x)   .  dU(x)  . 

aCx)  =  k(x)  — z where  the  value  of  — ; is 

dx  dx 

given  above. 

2 .   Linear  Displacement  Model  and  Linear  Properties 

In  the  first  example  the  exact  stress  function  was 

known  and  was  used  to  calculate  G.T(X)»   Consider  the  re- 

N 

suits  if  the  stress  field  is  approximated  by 


O(x)  =  k(x) 


dU(x) 
dx 


(52) 


u     dU(x)  .  .  ,    . 

where  — r-5 — -  xs  a  constant  over  each  element, 
dx 

This  is  the  constant  strain  approximation.   Since 

a  constant  strain  triangle  (CST)  finite  element  computer 

program  is  used  to  provide  the  O^, '  s  in  the  two-dimensional 

problems  considered  later,  using  -= —  =  constant  provides  a 

dx 

guide  to  the  error  which  may  be  expected. 

U(x)  as  given  in  Equation  (36)  is  calculated  at 
each  node  and  plotted  vs  x  below  assuming  linear  displace- 
ment between  nodes.   The  slope  of  each  line  segment  is  -; — 

r  dx 

and  is  shown  in  parentheses. 
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U(x) 


U(x)(|) 


Figure  11.   Assumed  Linear  Displacement  Model 


The  conjugate  nodal  components,  a  (x)  ,  are  now 


given  by 


ak 


ak    T 
°N(X)  B"f/   t(1+x)  dT]  ^N(X)  dx  (53) 

The  results  of  the  integration,  disregarding  the 


factor,  are 


1    4 
°1  =  6 


2    21 


3    50 


rr1      8 

°2    =  6 


2    24 
°2  =  T 


^3    55 
°2   =   T 


(54) 


The  conjugate  function  representation  of  stress  is 
given  by  Equation  (21)  and  the  results  are  plotted  on 
Figure  12  and  included  in  Table  1. 
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Figure  3,2.   Exact  Stress  Distribution  and  Approximations 
to  the  Exact  Stress. 
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TABLE  I.   Summary  of  Results  for  Problems  III  Al  and 
III  A2. 


0(X) 

x  Exact 

a(x) 
CBMl 

Error 

% 

a(x) 
CBM2 

Error 

% 

C.F.E. 

a(x) 
C.Ave 

Error 
% 

0.  0 

1.  4 

2.  12 

3.  24 

-0.33 
3.43 
11.7 
23.6 

14.25 
2.50 
1.66 

-0.13 

4.26 

12.06 

21.46 

6.5 

0.5 

10.58 

1 
2,6 
9,15 
20 

1 

4 

12 

20 

0 

0 

16.67 

NOTES:   1.   a(x)  CBMl  used  known   (x) 

2.  a(x)  CBM2  used  -z —  =  constant 

dx 

3.  cr(x)  C.F.E.  is  the  conventional  solution 
shown  in  Section  III  Alh 

4.  O(x)  C.Ave  is  the  average  at  a  node  of 
the  C.F.E  solutions 


5. 


Error  %  is  Exact-Approximate    -  nr. 

*  r   x  100 

Exact 
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Figure  (12)  shows  three  different  approximations  to  the 
exact  solution.   The  conjugate  basis  method  with  the  exact 
stress  distribution  gave  the  best  solution  at  node  4  and 
had  the  least  mean  square  error  throughout.   The  next  best 
approximation  was  the  conjugate  basis  method  where  the 
exact  stress  distribution  was  not  known.   The  conventional 
finite  element  solution  is  also  plotted. 

3 .   Linear  Displacement  Model  and  Constant  Properties 

This  is  the  same  problem  considered  in  the  previous 
section  only  this  time,  k(x)  =  k   =1. 


Equation  (53)  becomes 

ak     1 
~e /  \         °   fdU   .  e  .  N 

Vx)  ="T  J  dl  Vx)  d: 

The  nodal  stresses  are 


(56) 


111 
°1  =  °2    '      2 


2     2     3 
ai  =  °2  =   2 


3     3     5 
°1  =  °2    =   2 


(57) 


The  conjugate  basis  representation  of  the  stress 

field  is  found  using  Equation  (21)  and  is  shown  on  Figure 

13.   The  exact  solution  in  this  case  is  a(x)  =  2x.   The 

conventional  finite  element  solution  is  also  shown.   The 

ak 
factor  -  — r —  has  been  disregarded. 
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a  Exact 
x 


Figure  13.   cr(x)  vs  x  for  k(x)  =  1. 

4 .   Uniform  Strain 

Consider  the  same  three  element  models  with  k(x)  = 
600  to  cancel  the  factors  1/40  and  1/15  resulting  from  the 
numerical  calculations  and  U(x)  =  yr-  .   This  is  a  case 
where  node  4  is  displaced  0.3  units  along  the  x  axis. 
dU(x) 


dx 


=  constant  =  .1    In  this  case  the  exact  solution, 


the  conventional  finite  element  solution,  and  the  conjugate 
basis  solution  are  exactly  the  same.   This  is  shown  in  Fig- 
ure 14.   This  problem  corresponds  to  a  uniform  bar  under  a 
concentrated  end  load  P  =  30.  (EA) 
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X 


30 


Exact 

C.F.E. 

CBM 


Figure  14.   o(x)  vs  x  for  Uniform  Strain. 

5  .   Discontinuous  Displacement 

In  the  previous  examples,  the  stress  was  in  fact 
continuous  throughout  the  model.   The  CBM  gave  the  best  fit 
at  the  extremes  and  the  least  mean-square  error  over  all. 
However,  in  problems  where  the  stress  is  discontinuous  due 
to  a  discontinuous  displacement  approximation,  the  CBM  will 
not  be  the  best  fit.   However,  we  consider  examples  with 
stress  discontinuity  to  observe  how  the  CBM  tries  to  mini- 
mize the  total  error  in  a  least  squares  sense  in  various 
models.   The  direct  relation  between  discontinuous  load  and 
stress  obtained  here  is  true  for  the  in-plane  and  bar  prob- 
lems.  However,  other  problems  in  mechanics  may  have  other 
stress-load  relations.   For  example,  for  bending  structures 
(beams  or  plates),  discontinuous  stress  arises  from  discon- 
tinuous moment,  and  not  discontinuous  force. 


A3 


a.   Displacement  at  Node  3 

For  the  same  three-element  bar  given  a  displace- 
ment of  0.1  units  at  node  3,  the  displacement  function  is 


now 


U(x)  - 


x 

20 


0  <  x  <  0 


20    2  1  x  1  3 


(58) 


The  strain  for  elements  one  and  two  is  a  con- 
stant 0.05  and  for  element  three  is  zero.   The  exact  stress 
distribution  is  given  by  o(x)=  k(x)— -p — -  where  k(x)  =  600. 

The  conjugate  nodal  values  are  computed  using 
Equation  (23)  and  are  shown  below. 


~1    ~1    ~2    ~2    600 
°1  =  °2    =    °1    =  °2  =  TO 


~3    ~3 


(59) 


The  conjugate  stress  field  is  computed  using 
Equation  (21)  and  is  shown  on  Figure  15  along  with  the 
exact  and  conventional  solutions. 
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Figure  15.   o(x)  vs  x  for  Displacement  at  Node  3. 

The  conjugate  basis  method  spreads  the  error 
throughout  the  model  to  obtain  the  best  least  squares  fit, 
thereby  giving  a  stress  field  in  the  element  where  there  is 
no  stress.   Consider  the  area  under  the  stress  curve  as  a 
measure  of  the  average  error.   For  the  exact  solution  A  = 
30  x  2  =  60  units;  for  the  conjugate  basis  solution  A  = 
60.4  units.   This  gives  a  total  percentage  error  of  0.66%. 
However,  at  node  3,  the  CBM  value  is  only  53.3%  of  the  maxi- 
mum stress  at  that  point.   Clearly  the  conventional  finite 
element  method  gives  a  better  approximation, 
b.   Displacement  at  Node  2 

If  the  model  had  been  displaced  from  node  2 
instead  of  node  3,  the  results  would  be  similar  in  that  the 
conventional  constant  strain  element  would  be  a  better 
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approximation  and  the  CBM  spreads  the  error  over  the  entire 
model.   The  results  are  shown  below. 


40. 


30- 


Gx(PSI)  20- 


10- 


a   Exact 
x 


-10 


Figure  16.   o(x)  vs  x  for  Displacement  at  Node  2. 

Again  comparing  the  areas  under  the  stress 
curves  A(exact)  =  30;  A(CBM)  =  30.89.   The  error  is  2.96% 
overall  and  53.33%  less  than  the  maximum  at  node  2. 

B.   TWO-DIMENSIONAL  IN-PLANE  PROBLEMS 

The  development  of  the  computer  program  which  was  used 
to  solve  the  two-dimensional  in-plane  problems  is  discussed 
in  Section  IV.   Appendix  B  gives  a  description  of  the  dif- 
ferent finite  elements  which  were  used  to  obtain  results  to 
compare  with  the  results  obtained  by  the  CBM.   The  computer 
program  itself  is  listed  in  Appendix  C. 
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The  basic  model  used  to  verify  the  computer  program  was 
a  10  x  10  x  1  plane  divided  into  either  four  or  eight  ele- 
ments as  shown  below.   The  units  used  for  numerical  calcula- 
tions are  inches  and  pounds.   The  figure  also  shows  the 
global  node  numbers  and  the  element  numbering  system. 


10M 


4            / 
/        3 
/      6 

\         1 
2    \ 

5               \ 

\v       6 
5        \ 

7        / 
/            8 

Model  A 


Model  B 


Figure  17.   Typical  Two-dimensional  Models. 

Each  node  has  two  degrees  of  freedom  (DOF)  ,  which  are 
in-plane  displacements  u  and  v,  and  for  a  typical  node  n 
they  are  numbered  as  shown  below. 
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n 


D0F(2n)  =  v   =  v.    N 
n     (x,y) 


n 


D0F(2n-l)=  u   =  u 


n     (x,y) 


n 


Figure  18.   Reference  System. 

The  following  information  is  pertinent  to  all  of  the 
problems . 

E  =  30  x  106  lbs/in2 
y  =  0.3 
All  stresses  plotted  are  a  (x,y)  x  10   lbs/in 
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1 .   Model  A  with  Uniform  Edge  Displacement 


329 


329 


u(0,y)  =  0. 

v(0,y)  =  0. 

u(10.,y)  =  1. 

v(10.,y)  =  0. 


329 


Figure  19.   G  (PSI  x  10  )  for  Problem  III  Bl. 


The  plane  is  prevented  from  displacement  in  the  v 
direction;  therefore  c      =  0.   Equation  (A-l)  shows  that  the 
stress  will  be  uniform  across  the  model  since  the  strain  is 
uniform.   This  is  the  result  that  the  CST  solution  gives. 
Because  the  CST  solution  is  continuous  for  all  elements,  the 
conjugate  basis  solution  cannot  improve  upon  it.   Both  solu- 
tions in  this  case  are  exact,  for  plane  stress. 

This  problem  is  analogous  to  the  three  element  bar 
which  was  displaced  at  node  4.   See  Section  III  A4  .   In  the 
following  figure,  the  stresses  obtained  by  the  conventional 
CST  finite  element  solution  are  given  in  each  element,  while 
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the  stresses  obtained  by  the  conjugate  basis  method  (CBM) 
are  the  values  given  at  the  nodal  points. 

2 .   Model  B  with  Uniform  Edge  Displacement 

This  problem  is  analogous  to  both  the  problem  of 
the  bar  pulled  uniformly  at  its  end  (III  A4)  and  the  first 
Model  A  plane  problem  (III  Bl) .   It  confirms  that  mesh  A 
was  sufficient  for  the  problem  and  that  the  CST  is  a 
sophisticated  enough  element  for  this  type  of  problem.  The 
difference  in  boundary  conditions  between  this  problem  and 
problem  III  Bl  accounts  for  the  difference  in  stress  values 
In  this  case  c   <  0.  that  is  contraction  in  the  v  direction 

y 


is  allowed,  thereby  reducing  a 


x 


300 
Pi 


30Q 


300 


300     / 
/       300 

\^    300 

300   >v 

300        \ 

\300 
X    300 

300   \. 

i  ,, - 

300    / 
/         300 

300 


300 


u(0.,y)  =  0. 
v(10.  ,y)=  1. 


300 


300 


Figure  20.   a   (PSI  x  10  )  for  Problem  III  B2. 
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3 .   Model  B  with  Uniform  Edge  Displacement,  Modified 
Analytical  values  of  a   were  calculated  for  x  = 
x(10,y),  x  =  x(5,y)and  x  =  x(0,y)using  the  plane  stress 
equations  (A-l) .   The  values  of  strain  were  approximated  by 
dividing  the  total  contraction  in  the  v  direction  for  a 
given  y  by  the  original  length  10.   The  strain  e/  was  con- 
sidered constant  along  a  constant  y.   The  displacements 
were  calculated  in  the  CST  routine. 


The  values  of  O       calculated  were 

x 


a  (0.)  =  329 
x 


a  (5.)  =  292 
x 


a  (10.)  =  300 
x 


For  x  =  0  and  x  =  5  the  CBM  gives  the  better  approxi- 
mation.  At  x  =  10  the  difference  in  percentage  of  error 
between  the  CST  solution  and  the  CBM  is  less  than  0.67%. 
It  is  also  interesting  to  note  that  the  average  of  the  CST 
values  at  nodes  2  or  8  is  the  same  as  the  CBM  value.   This 
gives  mathematical  substance  to  the  engineering  practice 
of  using  the  average  value  at  a  node. 


u(0.,y) 
v(0. ,y) 
u(10.,y)  =  1 


51 


4 .   Model  A  with  Linear  Edge  Displacement 

Here  the  right  edge  is  given  a  linear  displacement 

as  shown  below.   Using  the  analytical  procedure  discussed 

4 
in  Problem  III  B3,  the  values  of  O       range  from  329  x  10 

x 

along  the  edge  (x,10)  to  zero  along  the  edge  (x,0).   The 
CBM  values  are  the  best  continuous  approximation  available 
for  this  model.   A  plot  of  a  (0)  vs  y  along  (0,y)  shows 
how  the  CBM  is  better  in  approximating  the  maximum  value. 


u(0. 

,y)   = 

0. 

v(0. 

,y)   = 

0. 

u(10 

.,y) 

=  y/io 

v(10 

-,y) 

=    0. 

24 


24 


Figure  22 


ax(PSI  x  10  )  for  Problem  III  B4. 
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Figure  23.   a   (0)  vs  y  for  Problem  III  B4 . 


A  comparison  of  the  percentage  of  error  at  (x,y)  = 
(0.,  10.)  that  the  CBM  value  is  7.4%  too  low  while  the  CST 
average  value  is  28.5%  low  and  the  minimum  CST  value  is 
50%  low. 

5 .   Model  B  with  Linear  Edge  Displacement 

The  applied  displacements  in  this  problem  make  it 
similar  to  Problem  III  B4;  however,  the  other  boundary  con- 
ditions in  this  problem  allow  for  movement  in  the  v  direc- 
tion.  In  part  the  displacement  solution  of  the  CST  program 
shows  that  at  x=0.,  x=5  .  ,  and  x=10.,  £   <  0.   Therefore, 
along  the  bottom  edge  where  ex  is  very  small,  the  effect  of 

E   becomes  dominant  and  a  negative  a   results, 
y  x 


53 


From  the  results  of  the  CST  displacement  solution, 

9u 
the  value  of  e   =  7: —  along  the  edge  (x,0.)  is  0.0026  and 

X       0  X 

9v 

the  value  of  e   =  tt—  ■  0.016  along  x=10.  and  0.015  along 
y    3y 

(5,y)  and  (0,y).   Using  these  values  of  £   and  the  average 
of  the  e     's  in  Equation  (A-l)  for  stress  yields 
a  (x,0)  =  -  73  x  103  PSI 

The  CST  solution  is  attempting  to  reflect  this  in 
elements  5  and  8.   The  CBM  does  show  the  compressive 
stresses  and  of  course  is  continuous  throughout  the  model. 
The  average  of  the  CST  values  at  nodes  where  more  than  one 
element  are  incident  is  close  to  the  CBM  value  at  that 
node.   For  example,  at  node  5  the  CST  average  is  150  and 
the  CBM  value  is  148.   This  again  confirms  the  existing 
practice  of  averaging. 


398 


152 


u(0. ,y)   =  0. 
u(10.  ,y)  =  y/10 


Figure  24.   o   (PSI  x  10  )  for  Problem  III  B5 . 
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6 .   Model  B  With  Piecewise  Linear  Edge  Displacement 

This  problem  has  boundary  conditions  similar  to 

problem  III  B  A;  but,  since  it  has  two  elements  along 

x  «  x(10,y)  and  the  applied  displacement  is  only  at  node  1, 

the  displacement  profile  of  the  edge  is  only  piecewise 

linear  with  discontinuous  slope. 

The  expected  stress  profile  is  greatest  at  (x,10) 

where  £   is  maximum  and  compressive  at  (x,0)  where  £  =0 
x  r  '  x 

but  a  compressive  £   exists.   Likewise  £   is  higher  along 

y  x       &         b 

(0,y)  because  £  =0,  than  along  (5,y). 

The  CBM  values  conform  with  the  expected  distribu- 
tion and  give  a  continuous  stress  field  which  is  the  best 
approximation  in  this  case. 


192 


u(0. ,y)  =  0 
v(0.,y)  =  0 
u(10.  ,0.)=  0 
v(10.  ,0.)=  0 
u(10. ,10.)=1 
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Figure  25.   a   (PSI  x  10   for  Problem  III  B6 


55 


7 .   Model  A  With  Center  Node  Displacement 

For  this  problem  the  stress  is  discontinuous  due  to 
the  displacement  of  node  5  in  the  interior  of  the  model. 
The  CST  solution  gives  an  exact  representation  of  the  stress 
field.   Element  1  is  in  compression  and  element  3  is  in 
tension. 

The  CBM  values  are  exact  at  the  nodes.   A  stress 
plane  constructed  linearly  thru  the  nodes  would  be  in  error 
at  all  'locations  other  than  the  nodes. 

If  a  three  dimensional  plot  is  made  with  stress  as  the 
vertical  axis,  the  volume  under  the  CST  plot  would  equal 

the  volume  under  the  CBM  plot  at  25  x  O 

max 
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-329 
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Figure  26.   a   (PSI  x  10  )  for  Problem  III  B7. 
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8 .   Model  A  With  Center  Node  Displacement,  Modified 

The  model  is  subject  to  a  discontinuous  stress  field 
due  to  the  internal  displacement  of  the  cneter  node.   If 
the  model  allowed  more  flexibility,  one  would  expect  the 
stress  to  be  discontinuous  in  the  region  of  node  5  and 
continuous  along  the  edges.   The  CBM  does  not  show  the 
abrupt  change  at  node  5  but  does  approximate  the  solution 
expected  of  a  more  flexible  and  hence  more  accurate  model. 
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Figure  27.   a   (PSI  x  10  )  for  Problem  III  B8 
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9 .   Model  B  With  Center  Node  Displacement 

This  problem  is  similar  to  the  previous  problem  on 
model  A  with  the  exception  that  the  boundary  conditions  on 
this  problem  allow  displacements  in  the  v  direction.   The 
stress  is  discontinuous  in  the  region  around  the  center 
node,  node  5.   However  it  is  unlikely  that  the  stress  is 
as  severely  discontinuous  as  the  CST  solution  indicates  in 
this  region.   Intuitively  a   has  a  maximum  tensile  value  at 
(0,5)  and  a  maximum  compressive  value  at  (10,5)  and  decreases 
in  absolute  value  away  from  the  line  of  action  of  the 
displacement . 

For  this  model,  subject  to  a  point  discontinuity, 
the  CBM  gives  a  better  overall  representation  of  the  stress 
field. 
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u(0.  ,y)  =  0 
u(5.  ,5.)=  1 


Figure  28.   c?x  (PSI  x  10  )  for  Problem  III  B9 
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10 .   Model  B  With  Uniform  Centerllne  Displacement 

This  is  another  problem  with  a  discontinuous  stress 
resulting  from  an  internal  displacement.   Note  that  to 
minimize  deflection  in  the  v  direction,  the  boundary  condi- 
tion v(5,5)  =  0  has  been  imposed. 

The  expected  value  of  a  (x)  for  x  >  5  is  zero.   The 
CST  solution  is  superior  from  two  standpoints.   It  reflects 
the  stress  discontinuity  very  well  and  it  gives  a  better 
approximation  to  the  stress  field  over  the  right  half  of 
the  model.   It  is  interesting  to  note  that  the  average  CST 
value  at  nodes  2  and  8  is  the  CBM  value. 
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Figure  29. 


a   (PSI  x  10  )  for  Problem  III  B10 

x  ' 
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11 .   Model  B  Modified,  With  Uniform  Centerline  Displacement 
This  mesh  is  a  modification  of  model  B  and  was 
designed  to  determine  how  rapidly  the  CBM  stress  would  con- 
verge to  zero.   Each  node  along  x=5  and  x=10  was  displaced 
an  equal  amount  to  force  rigid  body  motion. 

The  CST  model  responds  to  the  discontinuity  quickly 
as  expected.   The  CBM  does  try  to  converge  and  is  a  better 
approximation  over  this  part  of  the  model  than  the  approxi- 
mation with  the  previous  mesh  (III. B. 10).   The  conclusion 
is  that  if  a  discontinuity  occurs  across  an  entire  model, 
and  there  is  sufficient  room  for  the  continuous  function  to 
correct  itself,  and  the  area  of  interest  is  not  at  the 
discontinuity,  the  CBM  will  be  the  best  approximation  for 
a  continuous  field. 


187    -22     37 
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u(10.,y)  =  1 
u(5.,y)  =  1 
u(0.  ,y)   =  0 
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Figure  30.   a   (PSI  x  10  )  for  Problem  III  Bll 
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12 .   Model  B  Tapered  With  Uniform  Edge  Displacement 

This  model  is  another  modification  to  model  B.   The 
length  has  been  doubled  and  a  taper  introduced  by  making  the 
right  edge  one  half  the  height  of  the  left  edge.   Symmetry 
about  the  horizontal  centerline  has  been  maintained.   The 
model  simulates  a  tapered  plane  and  is  the  two  dimensional 
analog  of  the  one  dimensional  tapered  bar.   Both  the  CST 
solution  and  the  CBM  give  the  expected  stress  distribution. 
However  the  CBM  is  continuous  over  the  entire  model  and 
therefore  is  a  better  approximation. 
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u(20.,y)=  1 
u(0. ,y)  =  0 
v(0.,y)  =  0 
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Figure  31.   a   (PSI  x  10  )  for  Problem  III  B12 
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13 .   Composite  Material  Models  With  Uniform 
Edge  Displacement 

This  problem  consists  of  using  the  CBM  with  a  CST 
input  on  the  composite  material  models  used  by  Lin,  Salinas, 
and  Ito  [Ref.  6]  and  comparing  the  results  for  stress 
convergence  and  computer  time  usage.   The  models  used  in 
Ref.  6  are  shown  in  figures  (32)  and  (33).   The  PLISOP  model 
shown  in  figure  34  has  been  used  by  Salinas  subsequent  to 
publishing  Ref.  6.   The  stresses  resulting  from  using  the 
CST  and  LST  routines  on  figures  (32)  and  (33)  and  the  PLISOP 
routine  in  figure  (34)  (see  Appendix  B)  are  shown  on  figures 
(35,  (36),  and  (37)  along  with  the  CBM  values. 

In  figure  (37),  along  the  edge  (5.2,y),  where  there 
is  only  one  material  present,  the  CBM  yields  a  smooth  curve. 
However,  where  two  materials  are  present  and  there  is  a 
sharp  change  in  stress  due  to  the  different  Young's  Modulus 
of  the  two  materials,  the  CBM  tries  to  smooth  out  the  change 
in  a  manner  similar  to  the  discontinuous  problems  shown 
earlier.   As  a  result  several  CBM  model  values  in  the  region 
of  the  change  in  material  are  not  in  agreement  with  values 
obtained  from  the  other  solutions.   The  results  are  in  good 
agreement  away  from  this  area  of  change. 

An  analysis  of  the  amount  of  computer  time  each 
program  used  is  shown  in  Table  II. 


62 


CO 

H 

P-i 

Z 

/— V 

o 

w 

P4 

to 

S 

O 

M 

w 

« 

►J 

.-J 

P-c 

w 

o 

vD 

CN 

CM 

O 

r^» 

r^ 

• 

• 

** 

oo 

vD 

CO 

o 

14-1 

CD 
00 
CO 
CO 


CO 
CO    W 

O 


Pn 
O 
Q 

O 

CO 


CM 
CO 
CM 


o 

CO 
M 

►-) 
P-i 


o-  o 
o 
•J 

W   CM 


CO 


O 


Csl 


CO 
CO 


CO 

w 

O 


Pn 
o 

Q 

<r 

CM 


CM 


O 
CM 


m 

CO 


CO 
CM 


m 


oo 

CO 


m 
oo 


en 


CM 
CM 


«3- 
CM 


M  ^ 


CO 
H 


in 
oo 


CO 


CO 
00 


m 


t3 

0 

cd 
H 

CO     CO 

►J     0 
CO 

4-1    43 

•H    O 
M 

P-, 


to 
c 


O   tH 
CJ   -H 

<u  a 

CO    -H 
CO 

G 

•H 
0) 

6 

•H 
H 

J-i 

co 

4-> 

B 
o 
u 


CO 

£ 

CO 

« 

S3 

w 

CJ 

t> 

ps; 

PS 

H 

Pi! 

CO 

O 

rH 

s 

tn 

fn 

w 

PS 

s 

O 

H 

•J 

W 

o 

fe 

<3 

<: 

S 

Pn 

o 

H 

hJ 

H 

►H 

P4 

J=> 

O 

H 

.J 

£ 

W 

CJ 

H 

< 

pi 

> 

hJ 

pq 

£> 

H 

o 

55 

<: 

:=> 

Pm 

O 

fe 

M 

CJ 

CO 

CJ 

H 

CO 
CJ 

>4-l 

o 
o 

CO 
•H 

CO 
CX 

o 

CJ 


63 


Figure  32."  Finite  Element  Grid  System,  124  DOF 
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Figure    33.       Finite    Element    Grid    System    242,    DOF . 
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Figure  35*.   (x,0.)(PSI  x  105)  for  Problem  III  B13 
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Figure  36.   ax(x=y)(PSI  x  10  )  for  Problem  III  B13. 
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Figure  37 


0x(5.2,y) (PSI  x  105)  for  Problem  III  B13 
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IV.   THE  COMPUTER  PROGRAM  DEVELOPMENT 

In  developing  a  computer  program  to  determine  stresses 

utilizing  the  conjugate  basis  method,  four  major  requirements 

had  to  be  met: 

i.   The  inverse  of  the  global  fundamental  matrix  is  needed. 

This  requires  forming  the  element  fundamental  matrices 

and  merging  them  into  the  global  fundamental  matrix 
itself . 

ii.   An'  approximation  to  the  stress  field  in  the  element  is 
required . 

iii.   The  conjugate  nodal  stress  values  for  each  element  are 
formed . 

iv.   The  conjugate  basis  approximation  to  the  stress  field 
must  be  assembled. 

There  are  many  existing  programs  which  satisfy  the 
second  requirement.   In  this  thesis  a  program  utilizing  the 
CST  written  for  use  in  Ref.  6  has  been  used.   It  was  chosen 
because  it  was  available;  it  uses  a  simple  two-dimensional 
element  —  the  three  noded  triangle;  and  problem  solutions 
using  this  element  were  available  for  comparison  as  to  the 
computer  time  used  and  the  accuracy  of  results.   It  was  not 
an  objective  of  the  thesis  to  develop  an  optimum  conventional 
finite  element  solution.   Any  other  CST  program  or  a  more 
sophisticated  program  could  have  been  used.   The  advantage 
of  the  CBM  is  that  it  makes  the  stress  at  a  node  a  function 
of  the  entire  model  and  not  just  the  element. 

Therefore  the  development  of  the  CST  program  itself  will 
not  be  discussed  with  the  exception  of  a  brief  comment  on 
the  constitutive  relation  matrix,  C.   For  plane  stress 
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problems  where  Cf   =  T    =  T    =  o 

z     xz     yz 


C  = 


1  -  V' 


1 

y 
o 


y  o 
1  o 
o   .5(l-y)J 


where   {a}  =  [C]{e} 

For  the  CST,  each  displacement  field  is  approximated  by 
a  linear  field  over  each  element,  hence  the  local  basis 
functions  are  area  coordinates.   The  local  fundamental 
matrix,  which  is  the  inner  product  of  each  of  the  basis 
functions,  is  therefore  dependent  on  the  total  area  of  each 
element.   The  integration  can  be  found  in  Felippa  [Ref.  7] 
and  the  resulting  matrix  is 


A 
12 


2  11 
12  1 
112 


After  the  area  of  each  element  is  calculated  in  the 
normal  CST  routine,  the  element  fundamental  matrix  —  EFM  — 
is  formed.   It  is  immediately  merged  into  the  global  funda- 
mental matrix  —  GFM.   The  merging  routine  is  described  in 
section  III.A.l.d. 

After  the  entire  GFM  is  formed,  it  is  inverted  using 
subroutine  SYINV.   This  is  the  most  time  consuming  addition 
to  the  CST  program.   See  Table  II.   Notice  that  for  the  92 
element  model  SYINV  used  20  seconds  but  for  the  200  element 
model  the  time  was  133  seconds. 
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The  conjugate  basis  nodal  stresses  are  calculated  in 
two  steps  at  the  completion  of  the  standard  CST  program. 
The  assumption  is  made  that  stresses  will  only  be  calculate 
at  nodes.   Since  <(>..  c  1  at  a  node  (see  section  II. C.  and 
section  III. A. I.e.)  the  values  of  the  local  conjugate  basis 
function  at  a  node  are  just  elements  of  the  GFM  inverse, 
in  accordance  with  equation  (44).   This  assumption  is  not  a 
serious  limitation  because  the  finite  element  mesh  can  be 
arranged  to  have  a  node  at  any  point  of  interest  or  the 
stress  can  be  calculated  by  interpolating  between  neighboring 
nodes.   However  it  greatly  simplifies  the  computer 
calculations . 

Recall  that  the  conjugate  nodal  values  of  stress  for 
each  element  are  given  by  Eq.  (23)  as 


aN  =  /   a(x)  ^N(x)dx 


The  CST  value  of  a(x)  is  constant  for  each  element  and  is 

L 

ty    (x) dx 

0 
A  (e) 
is  found  in  Ref.  7  and  is  — r —  for  each  N.   Therefore  the 

A(e) 
conjugate  value  for  each  node  of  element  e  is  — - —  a(e). 


*;•*?■*<.) 


(60) 


The  conjugate  function  representation  of  stress  is  given 
by  Eq.  (21)  as 


e  „,N 


*<x)  =  I      °*    ^(x) 
e  =  l 
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In  the  computer  program  a(x)  is  calculated  one  node  at 

a  time.   The  value  of  a„  is  multiplied  by  the  element  of  the 

N 

GFM  inverse  in  the  rov;  corresponding  to  the  global  node  that 
the  local  node  N  is  incident  to  and  the  column  corresponding 
to  the  global  node  being  calculated.   This  multiplication 
is  carried  out  for  each  element  and  the  resulting  sum  is 
printed  out  as  the  nodal  stress. 

The  modification  of  an  existing  conventional  finite 
element  program  to  perform  conjugate  basis  calculation  of 
stress  requires  the.  addition  of  three  simple  algorithms, 
i.   An  algorithm  which  forms  the  global  fundamental 

matrix  (GFM).   This  algorithm  may  be  located  in  the 
element  stiffness  DO  loop,  where  the  element  funda- 
mental matrices  (EFM)  are  individually  formed  and 
whose  elements  are  inserted  in  the  appropriate  place 
in  the  global  fundamental  matrix,  C.p. 

ii.   An  algorithm  to  invert  the  global  fundamental  matrix 

AT 
C.r,  and  obtain  C    (see  equation  (12).   C.p  is  a  banded 

symmetric  matrix  and  inversion  algorithms  are  available. 

(In  this  work  E.  Wilson's  SYINV  subroutine  was  used.) 

This  algorithm  is  located  after  C»r  is  formed,  and  is 

the  most  time  consuming  aspect  of  the  conjugate  basis 

method . 

iii.   An  algorithm  forthe  calculation  of  nodal  stresses 

associated  with  equations  of  the  conjugate  basis 

method.   This  algorithm  uses  the  conventional  stresses 

of  the  finite  element  method  and  must  therefore  follow 

the  conventional  stress  calculation  routine. 


A  listing  of  a  CST  finite  element  routine  is  given  in 
Appendix  C  with  the  CBM  modifications  shown  enclosed  in 
boxes . 
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V.   CONCLUSIONS  AND  RECOMMENDATIONS 

As  a  result  of  this  study  several  conclusions  can  be 
made  concerning  the  use  of  conjugate  basis  functions  with 
in-plane  problems. 

Where  the  model  is  not  subject  to  interval  discontinu- 
ities, the  CBM  gives  the  most  accurate  approximation  to  the 
stress  over  the  entire  model  for  a  given  displacement. 

If  there  is  a  point  or  localized  discontinuity,  the  CBM 
will  be  a  poor  representation  of  the  stress  near  the  dis- 
turbance but  will  be  a  good  approximation  away  from  the 
disturbance.   Mesh  refinement  can  localize  the  disturbed 
area.   It  is  understood  that  points  of  stress  discontinuity 
are  usually  points  of  interest  and  that  mesh  refinements 
generally  increase  computer  time,  therefore  the  use  of  the 
CBM  is  limited. 

If  there  is  a  stress  discontinuity  across  the  entire 
model,  the  CBM  is  a  poor  choice  if  the  model  is  too  small 
or  too  coarse  to  allow  the  CBM  to  settle  out.   Refining  the 
mesh  may  be  worthwhile  if  a  continuous  approximation  to  the 
stress  is  wanted  in  regions  away  from  the  discontinuity. 
However  the  additional  computer  time  needed  places  a 
premium  on  this  marginal  gain. 

If  a  large  number  of  problems  are  run  on  the  same  model 
and  the  global  fundamental  matrix  can  be  inverted  once  and 
stored,  the  additional  computer  time  per  problem  for  the 
CBM  solution  is  minimal.   The  advantage  of  a  continuous 
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stress  solution  can  be  had  without  using  the  more  sophisti- 
cated elements  with  nodal  strain  degrees  of  freedom. 

Finally,  at  nodes  where  the  stress  should  be  continuous 
and  several  elements  were  incident  at  that  node,  the  average 
of  the  CST  stress  values  at  that  node  agreed  very  closely 
with  the  CBM  stress  value.   This  gives  mathematical  sub- 
stantiation to  an  engineering  practice. 

It  is  recommended  that 

1.  The  CBM  be  used  for  in  plane  problems  without  interval 
discontinuities  because  the  resulting  stress  field  will 
be  continuous,  which  the  actual  stress  field  is,  and  the 
stresses  will  be  consistent  with  the  displacements. 

2.  Further  studies  with  the  CBM  should  be  carried  out  on 
other  classes  of  problems.   For  example,  plate  problems 
or  beam  problems  might  be  better  approximated  with  CBM 
solutions . 

3.  Additional  studies  be  carried  out  to  determine  the  cost 
effectiveness  of  incorporating  and  using  the  conjugate 
basis  method  in  existing  conventional  finite  element 
programs . 

A.   An  investigation  be  undertaken  to  determine  if  the  CBM 
can  be  modified  to  yield  stress  discontinuity  when 
appropriate . 

5.   A  study  be  made  to  establish  the  relation  of  domain 
modeling  to  CBM  analysis. 
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APPENDIX  A 
THE  IN-PLANE  PROBLEM 

In  this  discussion  a  plane  is  defined  as  a  body  bounded 
by  two  surfaces  of  zero  curvature  whose  thickness,  h,  is 
mush  less  than  its  surface  dimensions.   For  in-plane 
problems  only  axial  loads  N(x,y)  are  considered.   There  are 
no  lateral  loads  or  moments.   Displacement  boundary  condi- 
tions  can  be  of  the  type   au  +  bv  =  c   where  a,  b,  and  c 
are  specified  coefficients;  and  u  and  v  are  displacements. 
In  plane  forces  may  also  be  specified  on  portions  of  the 
boundary  where  displacement  is  not  specified. 

N 


11  N 

/— *•   x; 


Figure  A-l.   The  In-Plane  Problem 

Since  k  is  small  in  comparison  with  the  plane  dimensions, 

C   <<  a  »CJ    and   t   »T    <<  T    •   This  suggests  the  use  of 
z      x   y         zx   zy      xy  ° 

plane  stress  conditions,  that  is,  a  >  T   »  T    are  considered 

*  *  *       z '   zx*   zy 

to  be  zero . 

The  basic  stress  strain  relations  become  the  plane  stress 
equations . 


77 


1  -p 


2  (Ex  +  y  V 


a      =  tt  (ye  +  e  ) 

y      !_y2         x       y 


E 


xy    1  +  p    xy 


(A-l) 


The  strain  displacement  relations  are 


e 


e  = 


in 

3x 

3v 

3y 

-ii  (  |a  +  |J  ) 

dx    dy 


xy 


1  (    3u    3v  . 

2  (  37  +  3^  } 


(A-2) 


Substituting  Eqns .  (A-2)  in  Eqns .  (A-l),  the  partial 
differential  equations  for  the-in  plane  problem  become 


1  -  P 


r  3  u       3v  -i 
2    dx      3y 


1  -  P 


r    3u    3  v  -i 
2       3x    dy 


xy 


E        r   _9_U      3_V   -1 

2(1  +  y)  L  3y    3x  J 


(A-3) 
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Stress  resultants  N.N.N    having  the  units  of  force 

x '   y    xy 


2 

per  unit  length  (//  in/in  )  are  defined  as 


h/2 

N        f  ad: 

x      *  x 

-h/2 


N 


xy 


h/2 

/ 
•h/2 


T    dz 
xy 


N   = 

y 


h/2 

/   a, 
-h/2 


dz 


(A-4) 


These  last  expressions  imply  internal  discontinuity  in 
stress  is  associated  with  a  discontinuity  in  applied  internal 
in  plane-force. 

For  plane  stress  the  equilibrium  equations  are 


[  F   =  0  : 


do           8t 
x    xy_ 

ay 


3x 


9N     3N 


8x 


*2.  = 


9y 


(A-5) 


7  F   =  0  : 

L  y 


do  3t 

jl  + xy. 


3y 


8x 


8N     3N 
-*  +  -^  =  0 


3y 


8x 


(A-6) 


The  finite  element  method  yields  a  system  of  algebraic 
equations  which  are  obtained  from  the  above  partial 
differential  equations  by  the  usual  method  as  given  by 
Zienkiewicz  [Ref.  8]. 
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APPENDIX  B 
THE  FINITE  ELEMENT  METHOD 

An  extensive  description  of  the  finite  element  method 
can  be  found  in  Zienkiewicz  [Ref.  8].   However  to  facilitate 
understanding  of  this  project,  certain  features  of  the 
elements  will  be  presented. 

In  all  finite  element  programs  certain  rules  must  be 
followed  in  selecting  the  assumed  displacement  field. 

1.  The  finite  element  displacement  field  must  accommo- 
date rigid  body  displacements. 

2.  At  common  boundaries  two  elements  remain  in 
continuous  contact  with  each  other.   Neither  gaps  nor  over- 
lappings  are  permitted. 

3.  Constant  straining  patterns  must  be  included  in 
higher  order  descriptions  of  displacement. 

CONSTANT  STRAIN  TRIANGLE  (CST) 

For  the  CST  the  u  and  v  in-plane  displacements  are 
approximated  by  a  linear  field  over  each  element.   Since 
strain  is  the  first  derivative  of  displacement,  this  gives 
a  constant  strain  element-. 
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The  area  coordinates  used  in  the  program  are  defined  as 


2(x2,y2) 


.(xj^y^ 


1^  =  A1/A 


L2  =  A2/A 


L3  =  A3/A 


3(xo»yq) 


Figure  B-l.   Area  Coordinates 

With  these  coordinates  the  linear  displacement  field 
over  the  element  is  simply 


u(x,y)  =  L1  n1      +   L2  u2   +   L3  u3 


v(x,y)  =  Lx  v1      +      L2  v2   +   L3  v3 


(B-l) 
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LINEAR  STRAIN  TRIANGLE  (LST) 


'   I 


Figure  B-2.   Linear  Strain  Triangle 

The  LST.  has  three  nodes  on  each  side  and  each  displace- 
ment field  is  taken  as  a  2nd  order  polynomial  over  the 
element  domain.   The  displacements  u  and  v  for  each  element 
are  as  follows 


u(x,y)  =  Nxu   +  N2u2  +  N3u3  +  N  u4  +  N  u5  +  N  ug 


v(x,y)  =  NlV;L  +  N2v2  +  N3v3  +  N^  +  N5v5  +  N^ 


(B-2) 


where 


Nx  =  2L1  -  1 


N4  =  4LlL2 


for  a  typical  corner  node 
for  a  typical  mid  side  node 


The  N   are  called  interpolation  functions  for  the 
assumed  displacement  field.   Differentiating  the  displacement 
fields  gives  linear  strain  over  each  element,  hence  the  name 
linear  strain  triangle  (LST). 
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PLISOP 

The  PLISOP  program  was  coded  by  Professor  Gilles  Cantin 
and  his  students  at  the  Naval  Postgraduate  School,  Monterey 
It  performs  a  plane  stress  or  a  plane  strain  analysis  by 
the  finite  element  method  using  numerically  integrated 
isoparametric  elements.   These  elements  have  four  sides, 
each  of  which  may  be  either  curved  or  straight.   The  same 
shape  functions  which  define  the  geometry  also  are  used  to 
describe  the  displacement.   The  strains  will  be  a  function 
of  the  number  of  nodes  per  side  of  the  element.   The 
comparisons  made  in  this  thesis  were  to  a  cubic  element. 
This  isoparametric  element  uses  a  cubic  displacement  field 
and  yields  parabolic  strain  (stress)  fields  over  each 
element . 


Figure  B-3.   12  Node  PLISOP  Element 

All  of  these  elements,  the  CST,  LST  and  PLISOP,  with 
displacements  being  the  only  degrees  of  freedom, yield 
discontinuous  stresses  along  element  interfaces.   The  CBM 
gives  continuous  stress  along  element  boundaries. 
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For  each  of  the  conventional  finite  elements,  the  order 
of  the  strain  approximation  is  one  less  than  the  order  of 
the  displacement  approximation.   However  for  the  CBM  the 
displacement  approximation  and  the  strain  approximation  are 
of  the  same  order. 
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APPENDIX  C 
COMPUTER  PROGRAM 
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